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Abstract 

Consider a database most of whose entries are marked but the 
precise fraction of marked entries is not known. What is known is that 
the fraction of marked entries is 1 — e, where e is a random variable 
that is uniformly distributed in the range (0, eo) .The problem is to 
try to select a marked item from the database in a single query. If 
the algorithm selects a marked item, it succeeds, else if it selects an 
unmarked item, it makes an error. 

How low can we make the probability of error? The best possible 
classical algorithm can lower the probability of error to O (eg) . The 
best known quantum algorithms for this problem could also only lower 
the probability of error to O (eg) . Using a recently invented quantum 
search technique, this paper gives an algorithm that reduces the prob- 
ability of error to O (eg) . The algorithm is asymptotically optimal. 

1 Introduction 

Classical search algorithms are robust. If we reduce the problem size, the 
algorithm has fewer items to search and the performance of the algorithm 
will almost always improve. Quantum search algorithms depend on delicate 
interference effects, any change in parameters leads to significantly different 



"Research was partly supported by NSA & ARO under contract DAAG55-98-C-0040. 



2 THE PROBLEM 



results. For example, if we reduce the number of states in the database by 
a factor of four, a quantum search algorithm that would have previously 
succeeded, will now fail with certainty. Quantum algorithms usually exhibit 
oscillatory behavior in their performance characteristics. As a result of this, 
if there is some quantity that we want to maximize, e.g. the probability of 
the system being in marked states, we will need very precise knowledge of 
the problem parameters. Given such knowledge, it is easy to fine-tune the 
algorithm so that it achieves a probability of success of unity. What happens 
if we do not have this knowledge? This can be a serious handicap for a 
quantum search algorithm. This paper describes a way around this problem. 

The original quantum search algorithm ^ considered the problem of find- 
ing a marked item in a large unsorted database with minimum queries to the 
database. For this type of problem it is usually enough to be able to obtain 
the correct answer with a constant probability since the procedure can be re- 
peated a logarithmic number of times to drive the probability exponentially 
close to unity. In that problem a logarithmic factor was not significant since 
the quantum search algorithm gave a square- root asymptotic improvement. 
However there are other important problems where the additional queries 
create a significant overhead and need to be minimized, e.g. when we are 
limited to a single query and have to find the answer with a probability 
approaching unity. One field in which this type of problem occurs is in pat- 
tern recognition and image analysis where each query requires a lot of signal 
processing and the consequences of making an error are catastrophic. 

2 The problem 

Consider the situation where a large fraction of the items in a database are 
marked, but the precise fraction of marked items is not known. The goal 
is to find a single marked item with as high a probability as possible in a 
single query to the database. For concreteness, say some unknown fraction 
(1 — e) of the items are marked, with e uniformly distributed in the range 
(0, eo) with equal probability. The search algorithm returns an item, if it is 
a marked item the algorithm is said to succeed, else if it is unmarked, the 
algorithm is said to fail. 

Classically the smallest error that can be obtained is O (eg). To achieve 
this, uniformly guess an input and evaluate the function. If the function out- 
puts 1, you're done, otherwise guess again and cross your fingers (because you 
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have no more evaluations left). One cannot do better than O {e^) classically. 
In this paper we show that the probability of failure for the new scheme is 
O (eg) , whereas that of the best (possible) classical scheme and that of the 
best known quantum schemes are both O (eg). 

Before considering the specific problem mentioned above, let us describe 
a general framework. Consider the following transformation 

URsU^RtU\s) (1) 
Rs = l-{1- exp (^1)) \s) (s| , Rt = I-il- exp (z0) \t) {t\ 

U is an arbitrary unitary transformation, Rt & Rg denote selective phase 
shifts of the respective state(s) by |. Note that if we were to change these 
phase shifts from | to tt, we would get one iteration of the amplitude ampli- 
fication algorithm 0, 0- 

The next section shows that if U drives the state vector from a source (s) 
to a target {t) state with a probability of (1 — e), i.e. ||f/is||^ = (1 — e), then 
the transformation ^ drives the state vector from the source to the same 
target state with a probability of (1 — e'^) . The deviation from the t state 
has hence fallen from e to e^. Note that this is different from the amplitude 
amplification framework where the amplitudes were getting amplified; over 
here it is more convenient to present the results in terms of the probabilities. 

The striking aspect of this result is that it holds for any kind of deviation 
from the t state. Unlike the standard search (or amplitude amplification) 
algorithm which would greatly overshoot the target state when e is small 
(Figure 1); the new algorithm will always move towards the target. As shown 
in Section 5, this feature of the framework can be used to develop algorithms 
that are more robust to variations in the problem parameters. 

3 Analysis 

We analyze the effect of the transformation URgWRtU when it is applied to 
the \s) state. As mentioned in the previous section, R^ & R^ denote selective 
phase shifts of the respective state(s) by | (t for target, s for source). We 
show that if ||f/ts||^ = (1 — e) , then 

\\{t\URsU^RtU\s)\\' = {1-e'). (2) 
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Figure 1: The original quantum search algorithm is very sensitive to the 
fraction of marked states. For example, when the fraction is about 3/4, the 
algorithm fails, the algorithm of this paper will always produce an improve- 
ment. 



In the rest of this section, the greek alphabet 6 will be used to denote 
|. Start with |s) and apply the operations U, Rs,U\ Rt & U. If we analyze 
the effect of the operations, one by one, just as in the original quantum 
search algorithm ^ , we find that it leads to the following superposition (this 
calculation is carried out in the appendix): 

U \s) (e^^ + {e'' - 1)') + \t) Uts {e'' - l) . 

To calculate the deviation of this superposition from \t), consider the 
amplitude of the above superposition in non-target states. The probability 
is given by the absolute square of the corresponding amplitude: 

{l-\\Utsf) ||(e^'+l|f/t.f (e^'-l)')f • 
Substituting \\Uts\f = {I — , the above quantity becomes: 
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Figure 2: By setting tlie six-state ancilla, b, to the superposition 
^ (|0) + |1) + |2) + |3) + |4) + |5) u^) where u = exp (-f ) , we 
get a I phase-shift of the states for which F{x) = 1 relative to those for 
which F {x) = 0. [TU] gives a technique for implementing this with qubits. 



4 Existing Algorithms 
4.1 Classical Algorithm 

The best classical algorithm is to select a random state and see if it is a t 
state (one query). If yes, return this state; if not, pick another random state 
and return that without any querying. Note that this algorithm requires 
a single query, not two. The probability of failure is equal to that of not 
getting a single t state in two random picks since if either of the two states 
is a marked state, the algorithm will succeed. This probability is equal to 
e^. Since e is uniformly distributed in the range (0,eo). The overall failure 

probability becomes^yi^^ = |eo 



4.2 Quantum Searching 

Boyer, Brassard, H0yer and Tapp |1] first described in detail an algorithm 
that succeeds with probability approaching 1, regardless of the number of 
solutions (it's a classical algorithm that uses quantum searching as a subrou- 
tine; of course, it can be made fully quantum.) The first quantum algorithm 
to be able to search in a single query with a success probability approaching 
1 was given by Mosca [S]. 

Mosca observed that the quantum counting algorithm of PJ (based on 
the original searching algorithm) produces a solution with probability con- 
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verging to 1/2. One easily converts this to an algorithm with probability of 
success converging to 1/4. Thus by using this algorithm as a sub-routine in 
another quantum search, we get success probability converging to 1. (This 
appears to be based on the observation in [4 that an algorithm that succeeds 
with probability exactly 1/4 can be amplified to one with success probability 
exactly 1 using only one quantum search iterate). In other words, the tech- 
nique Mosca uses is to take a search algorithm that succeeds with probability 
1/4 — X and then use one quantum search iteration to map it to an algo- 
rithm that succeeds with probability (1 — 12X^ — 16X^). Using this scheme, 
if the fraction of marked states of the database is 1 — e, one can easily obtain 
a marked state with a probability of error of 1 — — je^ by means of a 
single quantum query. The overall failure probability in this case becomes 

- 4^0 16^0 

A recent quantum search based algorithm for this problem is by Younes 
et al [Oj (actually it is somewhat unfair to compare it to the other algo- 
rithms since this was specifically designed to perform well when the success 
probability was close to |, not close to 1). This finds a solution with a prob- 
ability of (1 — cos 6') (^^^^^^^j^ + aV'e' ) ' "^here q = number of queries and 
6 = arccose ((59) from [Hj). When q = 1, the success probability becomes: 
(1 — e) (1 + 4e^) , hence the probability of error becomes e — 4e^ + 4e^. The 

overall failure probability becomes j^o j^^^ '^'^'^ = (l^o — fco + ^o) 



5 New algorithm 

As in the quantum search algorithm, encode the N items in the database in 
terms of log2 N qubits. The algorithm consists of applying the transformation 
WR^WRtW to |o). W is the Walsh-Hadamard Transformation and is the 
state with all qubits in the state. After this an observation is made which 
makes the system collapse into a basis state. 

In order to analyze the performance of this algorithm, note that the 
algorithm is merely the phase shift transformation URgU^RtU applied to 
\s) which has already been analyzed in section 3. U is the W-H transform 
(W) and the state s is the state (state with all qubits in the state), then 
= 1 ~ where e lies in the range (0, cq) . Therefore after applying the 
transformation WR^WRtW to |o) , the probability of being in a non-t state 
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Figure 3: Comparison of the failure probability of the ir/S phase shift algo- 
rithm with a classical algorithm, [H], and [S], when the fraction of unmarked 
states e, varies between & 0.2. 

becomes e^, i.e. the overall failure probability becomes de ~ i^o 

The performance of the algorithm is graphically illustrated in Figure 3. 

6 Extensions 

6.1 Multi-query searching 

In practice, a database search would use multiple queries. The technique 
discussed above extends neatly to the multi-query situation. As described in 
[7j, the multi-query algorithm based on the ^ phase-shift transformation is 

able to reduce the probability of error to after q queries to the database. 

^9 + 1 

A classical algorithm reduces the probability of error to Note that q = 1 
gives the same results as in section 5. 

6.2 Optimality 

Both the single query, as well as the multi-query algorithms are asymptoti- 
cally optimal in the limit of small eg. This is separately proved in 
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6.3 Quantum Control & Error Correction 

Connections to control and error correction might be evident. Let us say 
that we are trying to drive a system from an s to a t state/subspace. The 
transformation that we have available for this is U which drives it from s 
to t with a probability ||f/ts||^ of (1 — e) , i.e. the probability of error in 
this transformation is e. Then the composite transformation URsU^RtU will 
reduce the error to e'^. 

This technique is applicable whenever the transformations U, W ,Rs & Rt 
can be implemented. This will be the case when errors are systematic er- 
rors or slowly varying errors, e.g. due to environmental degradation of some 
component. This would not apply to errors that come about as a result of 
sudden disturbances from the environment. It is further assumed that the 
transformation U can be inverted with exactly the same error. Traditionally 
quantum error correction is carried out at the single qubit level where individ- 
ual errors are corrected, each error being corrected in a separate way. With 
the machinery of this paper, errors can be corrected without ever needing to 
identify the error syndrome. This is discussed in [7j and |H]. 
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7 Appendix 

We analyze the effect of the transformation (3) of section 3 when it is applied 
to the \s) state. As in section 3, 6 will denote |. We show that if ||f/ts||^ = 
(1 — e) , then 

{t\URsU RtU\s) =[l-t^) (3) 

starting state |s) 
after RtU 



after U 



RtU \s) = Uts {e'' - 1) \t) + U \s) 
U RtU\s) = \s) + Uts {e''-l)U \t) 



after Rg : 

RsUR,U \s) = (\s) + Uts {e'' - l) U \t) 



+ {e''-l) {\s) + Uts {e^'-l)\s){s\U 

= \s) {e^' + WUtsf {e'' - 1)') + Uts {e'' - l) U \t) 
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after the final U : 



URsURtU \s) = U \s) (e^' + \\Utsf {e'' - l)') + Uts {e'' - l) \t) 

Calculating the error To estimate the error, consider the amplitude of 
the above superposition in non-target states (this is due to the por- 
tion U \s) {e'^ + \\Utsf {e'^ ' portion Uts {e'^ - l) |^) 
is clearly in the target state). The magnitude of the error probability 
is given by the absolute square of this error amplitude: 



(1 - \\u^ 



ts\ 



')||(<:"'+iic..ir 



Assume ||t/ts||^ to be (1 — e) , the above quantity becomes: 



— e 
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